Boiling of nuclear liquid in core-collapse supernova explosions 
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We investigate the possibility of boiling instability of nuclear liquid in the inner 
core of the proto-neutron star formed in the core collapse of a type II supernova. 
We derive a simple criterion for boiling to occur. Using this criterion for one of best 
described equations of state of supernova matter, we find that boiling is quite possible 
under the conditions realized inside the proto-neutron star. We discuss consequences 
of this process such as the increase of heat transfer rate and pressure in the boiling 
region. We expect that taking this effect into account in the conventional neutrino- 
driven delayed-shock mechanism of type II supernova explosions can increase the 
explosion energy and reduce the mass of the neutron-star remnant. 
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I. INTRODUCTION 



Detection of electron antineutrinos from SN 1987A (see [l|, |2j, |3j) confirmed the previous 
theoretical ideas about neutrinos playing crucial role in the core collapse of type II super- 
novae. According to well established estimates, only about one percent of the gravitational 



binding energy (or (1.4 ±0.4) x 10 51 erg for SN 1987A; see {4]) is released in the thermal and 
kinetic energies of the expanding ejecta; the remaining part is carried away by different types 
of neutrinos (Sj , which are known to be trapped effectively within nuclear matter during the 
last stages of core collapse (see 6j and references therein). Therefore, initially it did not 
seem hard to explain the observed values of kinetic energy of the expanding ejecta using 
neutrino-nucleon interaction as an effective channel of energy transfer from the emitted neu- 
trinos to nuclear matter 7J, |8j. However, concrete realization of this mechanism of supernova 
explosions encountered with great difficulties. Thus, it was realized that the prompt shock 

liseconds, 



wave, generated during the core collapse, fades away at the time scale of several mi 
losing its energy on nuclear dissociation 9], [h| and on the radiation of z/z/-pairs 



As a 



result, in simulations, the shock starts to move inward, pushed down by the infalling matter, 
leading to the formation of a black hole rather than to explosion. 



Considerable efforts were made to remedy this situation. Bethe & Wilson 12] showed that 
the shock wave could be revived about 100 msec after the core bounce if neutrino luminosity 
were sufficiently enhanced. The way of such enhancement of neutrino luminosity was opened 
after neutrino diffusion inside the center of the proto-neutron star (PNS) was discovered 
and found to be convectively unstable (see 14j, ll5j and references therein). Nevertheless, 
numerical simulations of the neutrino-driven shock revival still yield contradictory results. 



While explosions with kinetic energy up to 1.72 x 10 51 erg are obtained in 16j (0.5 • 10 51 erg 
in 15], 10 50 erg in 17], 0.94 • 10 51 erg in [l8]), there is also a number of simulations where 
explosions are not observed at all [3, [2^] . 

This inconsistency between the observational data and the results of core-collapse nu- 
merical simulations motivated some researches to look for alternative mechanisms. The 



magnetorotational mechanism, proposed by Bisnovatyi-Kogan 21|, can produce explosions 
with energies up to 0.61 x 10 51 erg at the time scale of ~ 0.5 sec after the core bounce 



22 



231 ] . The acoustic mechanism, developed by Burrows et al. 24|, also produces explosion; 



although its energy is rather uncertain in their numerical simulations, it develops for a time 
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range of several hundred milliseconds. Thus, even with new mechanisms and effects taken 
into account, the resulting energy release is marginally short of the observed values. 

One can conclude that all mechanisms of core collapse require some additional engine to 
provide the observable explosions. Such a new engine is the subject of the present paper. 
By analogy with the work of Bethe et al. [13j , we consider an additional type of transport 
of nuclear matter in the core which is different from diffusion and convection. It occurs in 
the form of boiling, i.e., first order transition between different phases of nuclear matter. 
We start with the observation that supernova matter at the inner core can exist in several 
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26, 



23,123,1301. 

Among them are the following (in the 



phases and their mixtures 
order of increasing density): 

• the spherical nuclei phase 

• the phase consisting of elongated nuclei (often called pasta) 

• the slab-like nuclei 

• the phase with cylindrical holes 

• the phase with charged microscopic bubbles, or cheesed phase (we will use this last 
term below not to confuse it with macroscopic bubbles of which we will speak later) 

• the homogeneous supernova matter 

It is natural to expect that several phase transitions can occur during the evolution 
of nuclear matter during as well as after the core collapse. Numerical simulations of the 
nuclear-matter phase transitions in supernovae were usually aimed at determining the ther- 
modynamic properties at the pre-bounce stage of collapse, and were needed to understand 
the development of the prompt shock wave. Such transitions taking place in a rapidly chang- 
ing environment during collapse can be called "short-term" phase transitions. In our paper, 
we discuss what we call "long-term" phase transitions, which occur in nuclear matter after 
the bounce under the condition of relative mechanical and local thermodynamic equilibrium. 
The main goal of this paper is to examine the potential importance of such "long-term" 
phase transitions on the dynamics of supernova explosion. 

In the "normal" case of mechanical and local thermodynamic equilibrium, different phases 
occupy the corresponding radial shells of the PNS, and this spatial phase picture evolves 
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continuously and relatively slowly in time. However, if heating of the nuclear matter (which 
is mainly due to diffusion of neutrinos) is sufficiently strong and inhomogeneous, it can lead 
to a condition similar to that in an ordinary teakettle. Specifically, a bulk of particular phase 
can overheat and become unstable with respect to phase transition, and small bubbles of the 
lighter phase can spontaneously appear in its volume. The bubbles will grow and, due to the 
Archimedean force, rise upwards. Henceforth, such a process is called boiling by analogy. In 
this paper, we argue that this process can provide more efficient mechanism of heating the 
outer parts of the PNS (compared with the neutrino-diffusion and convection mechanisms) 
and generate additional pressure wave. 

The paper is organized as follows. In Sec. HIl we give the basic thermodynamic description 
of the co-existence of various phases in supernova matter after the bounce. In Sec. IIIIl we 
derive the criterion of boiling to occur, and, in Sec. IIV| we provide our numerical estimates 

n 

for the values involved in this criterion using the tabulated equation of state from [311]. In 
Sec. IVJ we consider a simple model of the boiling mechanism and derive numerical estimates 
characterizing the efficiency of this process. We summarize and discuss our results in Sec. IVII 



II. PHASE EQUILIBRIUM 

Just after the core bounce, the inward movement of the inner core significantly decreases, 
and the prompt shock wave moves outwards losing its energy mostly on nuclear dissociation 
and radiation of uu-pairs. The material within and around the PNS approaches a state 
of mechanical equilibrium (while the velocities of the convective motion are much smaller 
than the appropriate first cosmic velocities). In contrast with the rapid collapse during the 
pre-bounce stage of contraction, local thermal equilibrium is a good approximation for the 
post-bounce stage. 

Because of the radial density gradient, the PNS at this stage has onion-like structure. The 
density of its inner core is higher than the nuclear saturation density, and the homogeneous 
matter phase is thermodynamically preferred. In the outer layers of the PNS, the density 
is much less than the nuclear saturation density, and matter exists in the form of ordinary 
nuclei. A number of intermediate phases can exist between these two shells. Adjacent 
phases are separated by the surfaces of coexistence. Here we derive a simple condition of 
coexistence of phases applied to the situation under consideration. 
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We can consider a phase transition characterized by fixed pressure p, temperature T, 
baryon number B, electric charge C and lepton number L. The conservation laws of the last 
three quantities imply the existence of the corresponding chemical potentials: Hb, He and 
Hl- We can determine their values from the chemical potentials of the supernova matter 
components (we suppose that the supernova matter consists only of neutrons, protons, 
electrons and electron neutrinos): 

Since the number of particle species is greater than the number of independent charges, we 
have one more relation on the chemical potentials (the so-called beta-equilibrium condition): 

fl p + fi e = fi n + LL V . (2) 

For fixed values of p and T, our thermodynamic system tends to a state with a minimum 



value of the Gibbs free energy (see, e.g., [321 ]) 



<& = ^HiNi = n B B + ncC + n L L. (3) 

i 

Since we are interested only in electrically neutral phases, we have 

C = 0, (4) 

and the second term in (jHJ) vanishes. The third term, lilL, does not change during the phase 
transition because the electron neutrinos interact with nuclear matter very weakly. 
Finally, the condition of equilibrium between phases 1 and 2 is 

A*in (Poj Tq, li v q) = /i 2 „ (po, T , la v q) , (5) 

where the subscript "0" marks the values of the parameters right on the interface between 
the two phases. 



III. CRITERION OF BOILING 



Local thermal equilibrium described in the previous section is continuously disturbed by 
the process of diffusion of electron neutrinos away from the central region of the PNS. This 
causes inhomogeneous heating of the bulk of nuclear matter. If this heating is sufficiently 
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strong, it can lead to an overheat of a particular phase, which may result in its boiling - 
emergence of bubbles of the adjacent lighter phase in its bulk. These bubbles will then raise 
and grow, effectively transferring heat and lepton number. Boiling is a particular case of a 
non-equilibrium first order phase transition. In this section, we derive a necessary condition 
of boiling in the form of a relation on the parameters of the supernova matter. 

According to Eq. ((Sj), in the process of external heating, a heavier phase 1 becomes 
metastable with respect to its transition to a lighter phase 2 if the condition is reached such 
that 

£tin(Po + Sp, T + ST, fi u0 + 5fi u ) > fi 2 n(po + Sp, T + 5T, fi u0 + 5fi u ) , (6) 

where the positive increments of the thermodynamic variables correspond to their radial 
gradients. From this, we obtain the condition 



2n 



dT 



dT 



dT 

dr 



+ 



+ 



<9/i 2 „\ 



In 



d P Jt.u,, V dp /Ttfh> 
In the important particular case of hydrostatic equilibrium, we have 



d/J> u 
dr 

dp 
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(7) 



dp 
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-Pm9 



(8) 



where p m is the mean density, and g is the local acceleration of free fall. Substituting this 
into ([7]), we obtain 



9/i 



2d 



dT 



dT 



dT 
dr 



d/i 



2d 
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p m g 



(9) 



IV. NUMERICAL ESTIMATES 



It remains to check whether the condition of boiling derived in the previous section can 
be realized in the usual supernova core-collapse. For this purpose, we estimate the values of 
the partial derivatives in (Q using the numerical simulations of the EoS of supernova matter 



given m 



3l|. 
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Phase 


Ub, fcri 3 


MeV 


p, MeV • fm 3 


T, MeV 


Hu, MeV 


Nuclei 


0.02 


-2.230 


0.1807 


4.00 


81.613 




0.04 


-1.953 


0.4556 


4.88 


105.363 




0.05 


-1.837 


0.6111 


5.17 


114.092 




0.06 


-1.818 


0.7738 


5.44 


121.699 


Cheesed 


0.05 


-2.660 


0.5615 


5.07 


114.149 




0.07 


-2.189 


0.9086 


5.65 


127.934 




0.09 


-1.967 


1.2726 


6.08 


139.143 




0.10 


-2.093 


1.4432 


6.23 


144.213 


Homogeneous 


0.10 


-2.536 


1.4241 


6.18 


144.530 




0.11 


-1.202 


1.7203 


6.64 


147.952 




0.12 


0.154 


2.2267 


7.09 


169.312 




0.16 


7.525 


4.227 


8.84 


188.917 



TABLE I: The parameters from 



311 ] used in the estimate of the partial derivatives in 



The authors of [31] tabulate all essential parameters describing the supernova matter 
in three phases: nuclei, cheesed phase and homogeneous matter. In Table [U we show the 
parameters which are required in our analysis. Because the required values of derivatives 
are not listed in 3l), we try to obtain them from interpolation. Namely, we use the finite 



differences: 

Hn(p2,T 2 , Hu2) - // n (pi,Ti, n ul ) ~ -q^-(P2 - Pi) + -t^t(T 2 -Ti) + -q^{^u2 -/ii/l) , (10) 

neglecting the higher derivatives. For each phase, we solve a system of three equations 
for three unknown variables. The numerical values obtained in such a manner should be 
considered as estimates. They are presented in Table [III 

To check whether the conditions of boiling are realized, we use the conventional values of 
the other parameters. Numerical simulations (see, e.g., 19||) give the following estimates: 



g « 1.0 X 10 14 cm ■ sec 2 



dr 



jrri 

'10—20) MeV • km -1 , — <0. (11) 

dr 



Using this data together with Table [TTJ we obtain the following conditions for boiling: 
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Phase 




( dfl n \ 


I dfl n \ 

JvT 


Nuclei 


-1.353 


-2.626 


0.125 


Cheesed 


0.666 


5.182 


-0.201 


Homogeneous 


2.87 


1.289 


-0.0317 



TABLE II: The results for the partial derivatives. 



for the transition between the nuclear and cheesed phase (p m ~ 0.8 x 10 14 g-cm 



dT 

dr 



< -(0.4-0.8) MeV-knT 1 ; 



• for the transition between homogeneous matter and cheesed phase (p. 
10 14 g • cm -3 ), 



(12) 

,„ ^ 1.6 x 



dT 
dr 



> -(1.0-1.5) MeV-km" 1 



(13) 



The difference in the inequality signs in estimates (fl2"l) and (TL3]) is connected with the 
difference in the signs of the coefficients of the temperature gradients in ([7j) or (Q for the two 
phase transitions under consideration. Note that conditions ( fl~2l) and ( fTBl are overlapping 
and complementary, so it is very likely that one of them is satisfied. We, therefore, can expect 
that boiling of nuclear matter can take place inside the supernova core. However, this only 
demonstrates the possibility of principle, calling for additional thorough investigation of this 
issue. 



V. A MODEL 

In this section, we construct a simplified model of phase transition between a heavier 
phase 1 and a lighter phase 2, which, for definiteness, we take to be the cheesed phase and 
the phase of nuclei, respectively. 

We assume spherical symmetry of the PNS, so that the cheesed phase is located in a 
spherical shell with radial coordinate from R — H to R. The free-fall acceleration on the 
upper boundary of this region is equal to 

GM 14 / M\ /10km\ 2 „ , . 

J= _ = 1 .33xl0»(-)(— ) cm.sec- (14) 

where M is the total mass inside the sphere of radius R. 
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We can estimate H using the condition of hydrostatic equilibrium: 

pgH^Ap, (15) 

where Ap is the dimension of the pressure interval in which the cheesed phase exists. Ac- 
cording to Table [H this pressure interval is approximately equal to 0.83 MeV-fm -3 . We 
thus have 

TT Ap „ „ / Ap \ /10 14 g-cm" 3 \ /1.33 x 10 14 cm • sec" 2 \ , 

H ps ps 1.0 X 5- km. 16 

pg V0-83MeV-fm- 3 y \ p J \ g J 

If we assume that boiling takes place in the whole volume of phase 1, then the total mass 
of the boiling matter can be estimated as 1 

M > - **P#H S 6.1 x 10- ( . 83M ^, fm -3 ) M . (17) 
The densities of phases 1 and 2 are related by p\ and p 2 = pi(l — e), where e <C 1. 



According to 34|, e is equal to 0.1 for the phase transition between the cheesed phase and 
homogeneous matter phase, and to 0.2 for the transition between nuclei and cheesed phase. 
Lower estimates for e are present in (m): it equals to 0.07—0.08 for both phase transitions if 
the entropy per baryon (which is conserved in the simulations) S/A = 1.0; for S/A — 1.5, the 
value e = 0.05 is obtained for the transition between nuclei and cheesed phase, and e = 0.01 
for the transition between cheesed phase and homogeneous matter. To be conservative, we 
use the lowest estimate in this paper, namely, e = 0.01. 

The maximum acceleration which can be reached by a raising bubble (neglecting the 
liquid resistance) is 

eGM 

«max — ■ (1°) 

The maximum velocity that can be reached by the bubble is then 



1 Recent simulations 33jj show that the total mass of different phases between nuclei and homogeneous 
nuclear matter (collectively called "pasta phases" in [331 ]) inside the supernova just before the bounce can 
amount to 0.13—0.30 M Q . This estimate somewhat differs from (|17|) because, before the bounce, one deals 
with small pressure gradients in an (almost) free-falling matter, while, after the bounce, we have matter 
close to hydrostatic equilibrium with its large pressure gradients. It is possible, in principle, that the 
pre-bounce boiling, if it occurs, can affect the propagation of prompt shock wave, but we do not consider 
this issue here. 
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Matter in the boiling volume will convectively move in both directions. We can expect 
that roughly half of matter moves upwards (the bubbles and the surrounding matter) and 
the other half moves downwards with the same average velocity (from the momentum con- 
servation). Therefore, the "connective boiling" overturn will be established. We can try to 
estimate the maximum efficiency of this overturn. 

The first effect to be discussed is heat transfer. If the bubbles fill half of the boiling 
volume, the heating rate at the surface of the boiling layer is 

Qmax ~ ^ x A-nR 2 v m ^ q 

= 7.7 x i(P ( JL_) 2 ( ) ( q ) erg . sec -i . (20) 

VlOkmy V 5100 km- sec" 1 / V 0.015 MeV-fm" 3 / 6 v ; 

Here, q is the specific volume heat of evaporation, and its value is estimated as q = 
0.015 MeV ■ fm -3 from [2_5J. The value (|20[) is comparable to the neutrino luminosity. This 
is, in principle, the upper estimate of heating rate which corresponds to maximal workload 
of the "boiling machine." In reality, the boiling heat transfer should work together with dif- 
fusion and/or convection below and above the boiling shell, increasing the net heat transfer 
rate. 2 This should lead to an increase in pressure behind the shock, providing more efficient 
conditions for shock revival. 

The second effect is the momentum transfer. The maximum mechanical pressure that 
the bubbles can exert is estimated as 

Pmax ~ pv 2 max = 2eAp = 1.7 x 10- 2 MeV • fm" 3 , (21) 

where the numerical value corresponds to the cheesed phase. This is much smaller than Ap 
since e <C 1; therefore, including this contribution to pressure does not seriously affect our 
previous calculation. But it can provide an additional barrier for the infalling matter to 
reach the inner core, reducing the final mass of the neutron star. 

We should admit that most of the above estimates represent upper limits corresponding 
to the maximal workload of the "boiling machine." In a subsequent paper, we will discuss 
these processes in greater detail. 



2 This can be compared to an electric circuit with a series of resistances. As one of the resistances in the 
series is shortened out, the total current increases. 



11 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we have demonstrated the possibility of boiling of nuclear liquid in the 
supernova core after the bounce. If it occurs, it can lead to effects increasing the efficiency 
of the neutrino-driven mechanism of supernova explosions. Among these effects are the 
following: 

• The increase of heat transfer rate from the inner core to the neutrinosphere. This 
increases the mean neutrino energy, making the delayed-shock mechanism more effi- 
cient. 

• The increase of pressure in the boiling region. It provides an additional barrier between 
the infalling matter and the inner core, thereby reducing the mass of the neutron-star 
remnant. 

We expect that taking into account the effect of boiling in the conventional delayed- 
shock/acoustic mechanism is important and will enable one to explain the energetics of the 
supernova explosions in a more simple and self-consistent way. 
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